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We apply a hard-core version of the Faddeev differential equations to the 3 He 4 He2 three-atomic 
system. Using these equations we calculate the binding energy of the 3 He 4 He2 trimer with the 
LM2M2 potential by Aziz and Slaman and more recent TTY potential by Tang, Toennies and Yiu. 

1 Introduction 

There is a great number of experimental and theoretical studies of the 4 He three- 
atomic system (see, e. g., GJ-G3 1 and references cited therein). The non-symmetric system 
3 He 4 He2 found comparatively little attention. We can only mention the recent works |6|, 
1 8 1, and 1 13 1 where the 3 He 4 He2 trimers were treated alongside with small 4 He clusters. 
Until now only the bound states of the 3 He 4 He2 system have been studied numerically. 
There are still no scattering calculations reported for this system. 

The 4 He trimer is known in particular for the Efimov's nature of its excited state (see 
l3l l6l l8l fl4l . The binding energy of the 4 He dimer is extremely small (about 1 mK) on 
the molecular scale. The large spatial extension of the 4 He2 bound state generates a long- 
range effective interaction between a 4 He dimer and additional 4 He atom which results in 
a possibility of existence of extremely extended 4 He three-atomic states. 

Being a more light particle than 4 He, the 3 He atom supports no bound state with the 
4 He counterpart and no 3 He dimer exists. Thus, the 3 He 4 He2 is even a more loosely bound 
system than the 4 He trimer. According to the hyperspherical adiabatic calculations of |6||8) 
and Monte-Carlo investigation of 1 13 1 the realistic He-He potentials such as LM2M2 1 15 1 
and TTY 1 16 1 support only one bound state of the 3 He 4 He2 trimer with the energy of the 
order of 10-15 mK. 

Notice that the 4 He/ 3 He three-atomic systems belong to the three-body systems whose 
theoretical treatment is quite difficult. The difficulty is mainly due to the two reasons. First, 
the low energy of the practically on-threshold bound states makes it necessary to consider 
very large domains in configuration space with a size of hundreds of A. Second, the strong 
repulsive part of the He-He interaction at short distances produces large numerical errors. 
Like |9 11 1, the present work is based on a mathematically rigorous hard-core version of 
the Faddeev differential equations. This method allows to overcome the strong-repulsion 
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problem. The first of the problems just mentioned is tackled by choosing sufficiently large 
grids. 

This note represents rather a first step in an extension of the numerical approach of 
(9] II II to the case of three-body systems including particles with different masses. In the 
nearest future we plan not only to continue our study of the 3 He 4 He2 bound state but 
also to perform calculations of the scattering of a 3 He atom off a 4 He2 dimer. Here we 
only outline the method employed and report our first results for the binding energy of the 
3 He 4 He2 system. 



2 Formalism 

In describing the 3 He 4 He2 three-atomic system we use the reduced Jacobi coordinates 
1171 x a ,y a , a = 1,2,3, expressed in terms of the position vectors of the atoms r; e R 3 and 
their masses m,-, 
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where (a, p , y) stand for a cyclic permutation of the atom numbers (1,2,3). The coordinates 
x a>ya fix the six-dimensional vector X = (x a ,y a ) £ R 6 . The vectors xp,yp corresponding 
to the same point X as the pair x a , y a are obtained using the transformations 



xp = cp a x a + sp a y a , yp = -sp a x a + cp a y a , 



(2) 



where 
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In the following we assume that the 4 He atoms are assigned the numbers 1 and 2 while the 
4 He atom has the number 3. By c we denote the hard-core radius which will be taken the 
same (in coordinates x a ) for all three inter-atomic interaction potentials. A recent detail 
description of the Faddeev differential equations in the hard-core model which we employ 
can be found in |9|. Nevertheless we outline here some essential characteristics of the 
hard-core Faddeev approach needed for understanding our numerical procedure. 

Since the 4 He atoms are identical bosons the corresponding Faddeev component F3 (X3 , y3 ) 
is invariant under the permutations the particles 1 and 2 which implies 



^3(-x 3 ,y3) = F 3 (x ?n y 3 ). 



(3) 
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The identity of the two 4 He atoms also implies that the Faddeev components F\ (xi ,yi) and 
^2(x2,y2) are obtained from each other by a simple rotation of the coordinate space. Thus, 
we only have two independent Faddeev components, the one associated with the 4 He- 4 He 
subsystem, Fi(x,y), and another one, say Fi(x,y), associated with a pair of 3 He and 4 He 
atoms. The resulting hard-core Faddeev equations read 



>P (a) (Xa,ya) 



0, (5) 
a= 1,3, 



where *FW and 'J'' 3 ' denote the total wave function ^{X) of the 3 He 4 He2-system written 
via the Faddeev components F\ and F3 in different coordinates xi,yi and X3,y3. More 
precisely 

y W (x hyi ) = Fi( Xl , yi ) 

+Fi(c 2 lXl+52iyi,-*2lXl +C 2 iyi) 
+F 3 (C31X1 +^3iyi,-i3ixi +c 3 iyi) 

and 

y {3 \x 3l y 3 ) = F 3 (x 3 ,y 3 ) 

+Fi (C13X3 + S13Y3, -513x3 + ci 3 y 3 ) 

+F\ (C23X3 + S23Y3 , -523X3 + C23Y3 ) • 

By Vi and V3 we denote the same interatomic He-He potential recalculated in the corre- 
sponding reduced Jacobi coordinates xi and X3. 

In the present investigation we apply the above formalism to the 3 He 4 He2 three-atomic 
system with total angular momentum L = 0. Expanding the functions F\ and F3 in a series 
of bispherical harmonics we have 

fix(x,y)^E y«o(5y), 1,3, (6) 

1 xy 

where x = |x|, y = |y|, x — x/x, and y = y/y. (Notice that by only the terms fj (x,y) 
with even momenta I are nonzero.) As a result the equations and boundary conditions 
are transformed to the following partial integro-differential equations 

(-£-£ + '<' +1 >(? + 7H^> 

0, x < c 

-V a (x)^ a Hx,y), x>c ' a=1 ' 3 ' t?) 
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and partial boundary conditions 

= 0, a=l,3. (8) 



y\ a) (x,y) 



fee) 

The partial wave functions , a = 1 , 3, read as follows 
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where (cf. |I7 ) 
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k max =-(l + X + l' + X'). 

Here Pk{f\) is the Legendre polynomial of order k. In the above, the standard notation for 
the 3- j, 6- j, and 9- j Wigner symbols, as defined in 111 81 . is used. We also use the notation 



*p<xCn) = \J Cp a -f 2 +2cp a sp a xyr| + s2 a 3;, 
ypa(ll) = \J s \ a x ~ 2c Pa s Pa J 0'1l + c j a y- 
We conclude the section with the asymptotic boundary condition for a 4 He 3 He 2 bound 
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as p = \Jx 2 +y 2 — ► °° and/or y — > «>. Here we use the fact that the helium dimer 4 He2 has 
a bound state and this state only exists for / = 0; £ c i stands for the 4 He2 dimer energy while 
\|/rf(x) denotes the 4 He2 dimer wave function which is assumed to be zero within the core, 
that is, \\fd(x) = forx < c. 



We employed the Faddeev equations Q, the hard-core boundary condition l|8}, and the 
asymptotic condition (II 01 to calculate the binding energy of the helium trimer 3 He 4 He2. 
As He-He interaction we used the semi-empirical LM2M2 potential of Aziz and Slaman 
1151 and the latest theoretical potential TT Y of Tang, Toennies and Yiu 1 1 6 1 . In our present 
calculations we used the value h 2 /m — 12.1 192 KA 2 where m stands for the mass of a 
4 He atom. (Notice the difference between this more precise value and the value h 2 /m = 
12.12 KA 2 which was used in the previous calculations IU031 ■) Both the LM2M2 and TTY 
potentials produce a weakly bound state for the 4 He dimer. We found that the 4 He-dimer 
energy is 1.309 mK in case the LM2M2 interaction and 1.316mK for the TTY potential. 
Both LM2M2 and TTY support no bound state for the 4 He 3 He two-atomic system. 

As in (5] II II we considered a finite-difference approximation of the boundary-value 
problem (Q|8]^J in the polar coordinates p = y/ 'x 2 + y 2 and = arctan(y/jc). The grids 
were chosen such that the points of intersection of the arcs p = p,, i = 1,2, . . . ,iVp and the 
rays = 8/, j = 1,2,... ,Nq with the core boundary x = c constitute the knots. The value 
of the core radius was chosen to be c = 1 A by the same argument as in II II . Also the 
method for choosing the grid radii p, (and, thus, the grid hyperangles 9 ; ) was the same as 
described in | Qj. 

In the present investigation we were restricted to considering only the two lowest par- 
tial components fQ l \x,y) and fQ 3 \x,y) and therefore we only dealt with the two partial 
equations of the system corresponding to / = 0. We solved the block three-diagonal 
algebraic system, arising as a result of the finite-difference approximation of (1711811101. on 
the basis of the matrix sweep method 1 19 1. This method makes it possible to avoid using 
disk storage for the matrix during the computation. 

The best possible dimensions of the grids which we employed in this investigation 
were N p = 600 and Nq = 605. We found that on the 600 x 605 grid with p max = 200 A the 
LM2M2 potential supports the bound state of the 3 He 4 He2 with the energy E t = 7.33 mK 
while the corresponding binding energy produced by the TTY potential is E t = 7.28 mK. 

Our figures for E t correspond to the lowest possible dimension of the system Q. 
We consider this as reason why our results show a significant underboundedness of the 
3 He 4 He2 trimer as compared to the available results for E t obtained for the TTY potential 
on the basis of the variational VMC (9.585 mK 03)) and DMC (14.165mK D3) meth- 
ods and for the LM2M2 potential on the basis of a one-channel hyperespherical adibatic 
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approximation of the Faddeev differential equations (9.682mK |8|) and (10.22 mK |6|). 
We think the situation will change when more partial waves in will be employed. A 
certain (but rather small) deepening of the binding energy E t may also be expected due to 
choosing the grids with larger A^e and A^ p . 
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